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We study the evolution of the equation of state of a strongly interacting quark system as a 
function of the diquark interaction strength. We show that for the system to avoid collapsing 
into a pressureless Boson gas at sufficiently strong diquark coupling strength, the diquark-diquark 
repulsion has to be self-consistently taken into account. In particular, we find that the tendency 
at zero temperature of the strongly interacting diquark gas to condense into the system ground 
state is compensated by the repulsion between diquarks if the diquark-diquark coupling constant 
is higher than a critical value Ac = 7.65. Considering such diquark-diquark repulsion, a positive 
pressure with no significant variation along the whole strongly interacting region is obtained. A 
consequence of the diquark-diquark repulsion is that the system maintains its BCS character in the 
whole strongly interacting region. 

PACS numbers: 12.38.Mh, 03.75.Nt, 26.60.Kp, 24.85. +p 



I. INTRODUCTION 

Shortly after the discovery of asymptotic freedom in QCD [l[, it was noted Q that the superdense matter might 
consist of weakly interacting quarks rather than of hadrons. Asymptotic freedom implies that at very high baryon 
density, QCD is amenable to perturbative techniques Q , and since the cores of neutron stars are formed by superdense 
matter, possible applications of those results to astrophysics Q were envisioned. Later on, however, it was understood 
that the ground state of the superdense quark system, a Fermi liquid of weakly interacting quarks, is unstable with 
respect to the formation of diquark condensates Q, a non-perturbative phenomenon essentially equivalent to the 
Cooper instability of BCS superconductivity. In QCD, one gluon exchange between two quarks is attractive in 
the color-antitriplet channel. Thus, at sufficiently high density and sufficiently small temperature T, quarks should 
condense into Cooper pairs, which are color antitriplets. These color condensates break the SU(3) color gauge 
symmetry of the ground state producing a color superconductor. In the late 90's the interest in color superconductivity 
(CS) was regained once it was shown, on the basis of different effective theories for low energy QCD [5j, that a 
color-breaking diquark condensate of much larger magnitude than originally thought may exist already at relatively 
moderate densities (of the order of a few times the nuclear matter density) and therefore it might be realized in 
compact stars. At densities much higher than the masses of the u, d, and s quarks, one can assume the three quarks 
as massless. In this asymptotic region the favored state results to be the so-called color-flavor-locking (CFL) state 
Q, characterized by a spin-zero diquark condensate antisymmetric in both color and flavor. 

Nevertheless, this picture breaks down at intermediate densities due to the mismatch between the Fermi momenta of 
different quarks produced by the strange quark mass M s and the constraints imposed by electric and color neutralities 
Q. That is, although the validity of the CFL phase at asymptotically large densities is well established, the next 
phase down in density is still a puzzle since as a consequence of the pairings with mismatched Fermi surfaces the 
phase exhibits chromomagnetic instabilities 0] ■ 

One possible scenario where this instability can be avoided occurs if in the region of moderate-low densities the 
strong coupling constant becomes sufficiently high (Gd ~ Gs ~ l/ 7 ^- 2 , with Gs and Gs denoting the diquark and 
quark-antiquark coupling constants respectively) @, Q. On the other hand, the increase of the coupling constant 
strength at low density can modify the properties of the ground state as indicated by the significant decrease of the 
Cooper-pair coherence l eng th, which can reach values of the order of the inter-quark spacing [1 01 ] . As already found in 
other physical contexts [ill, this fact strongly suggests the possibility of a crossover from a color-superconducting BCS 
dynamics to a BEC one |l2l|-[l5j. where although the symmetry breaking order parameter (the diquark condensate) 
is the same, the quasiparticle spectra in the two regions are completely different. As we will show by numerical 
calculations, in the BCS region, where the diquark coupling is relatively weak, the energy spectrum of the excitations 
has a fermionic nature, while in the strong-coupling region, formed by the BEC molecules, the energy spectrum of 
the quasiparticles is bosonic. 

As mentioned above, the combination of high densities and relatively low temperatures could exist in the dense cores 
of compact stars. The cores of neutron star remnants from supernovae collapse have densities several times larger 
than the saturation density of nuclear matter and temperatures several orders smaller than the superconducting gap. 
Under these conditions diquark pairs can form and resist the evaporation due to thermal effects. Then, it is natural 
to ask if a BEC of diquark pairs can take place at the moderately high density that the inner core of neutron stars 
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can reach. 

In this paper, we will show, through the mean-field analysis of the equation of state (EoS) of a simple system 
with four-fermion interactions, that by increasing the diquark coupling strength the matter pressure decreases up to 
negative values once the crossover from the BCS region to the BEC one takes place. This result would hint that if the 
density decreases to values where the coupling becomes sufficiently strong, the matter pressure turns to be negative 
and the system becomes unstable under the effect of gravity. Nevertheless, this is a naive picture that ignores the 
diquark-diquark interactions. As it was pointed out in Ref. [16j, together with the fact that there exists an attractive 
channel between quarks that favors the diquark formation, there is, as a corollary, a diquark-diquark repulsion. This 
repulsion is due to the cross-channel unfavorable correlations between the quarks belonging to different diquarks. 
Hence, when the diquark repulsion is self-consistently taken into account in the EoS of this system, the instability 
previously found in the strong coupling region is removed, and the pressure is stabilized with no significant variation 
through out that region. The increase of the diquark repulsion, which is produced by the raise of the energy gap in 
the strong coupling region, compensates the effect of the decay of the chemical potential, that as known, makes an 
important contribution to the EoS [13] ■ Yet, the price of the stabilizing effect produced by the diquark repulsion is 
that the Bose-Einstein nature of the strongly interacting system under study is lost, as we will discuss below. Our 
finding is calling attention on one hand to the necessity of including the diquark-diquark repulsive potential in the 
studies of the BCS-BEC crossover of strongly interactions, something that has been ignored up to now in previous 
works, and on the other hand, it is indicating that to include the diquark-diquark repulsion can completely change 
the understanding of this phenomenon in the context of strong interactions. 



II. BCS-BEC CROSSOVER AND QUASIPARTICLE SPECTRUM 

Our main goal in this section is to determine through a numerical analysis the threshold value of the attractive 
coupling constant between quarks that marks the crossover from the BCS to the BEC region in the frame of the 
simple model under consideration. 

For our investigation, we consider a simplified pure fermion system represented by the four-fermion interaction 
Lagrangian density of Ref. [l2[ 

C = Qiyfdn + 7oM - m)V> + |(^ 75 ^ T )(V' T Ci7 5 V') (1) 

In (fTJ, C = i"/oj2 is the charge conjugation matrix, m the fermion mass, fi the chemical potential defining the 
Fermi energy, and g the attractive coupling constant in the J p = + channel that parameterizes the strength of 
the interaction. Varying the strength of g yields the crossover from BCS (for a weak g) to BEC (for a strong g). 
The results we will obtain should not qualitatively change when additional internal fermion degrees of freedom, other 
than spin, are considered. It is due to the fact that the essence of the phenomenon under investigation is uniquely 
related to the change in the nature of the spectrum of the quasiparticles, which is determined by the variation of the 
diquark-pair binding energy as a function of the coupling constant strength. 

After introducing the Hubbard-Stratonovich transformation with gap parameter A = (gip T Cij5ip/2) , we have that 
the system free energy at finite temperature in the mean-field approximation is 

1 ~ r ei 3 k A 2 
^ = -7jE o 7 Trln t/ 3G "(-- k )] + y, (2) 

rp 

where G~ 1 (iu> n , k) is the inverse propagator in Nambu-Gor'kov space in the field basis ^ T = (ip, ipc), with ipe = Ctp 
being the charge-conjugate spinors, 

G^ 1 (iuj n , k) = (w„ + /icr 3 )7o - 7 • k - m + ij 5 Aa + + i7 5 A*cr_ (3) 

Here, w„ = (2n+l)7r//3 are the fermion Matsubara frequencies, and a± = ai±ia2, with 01,2 denoting the corresponding 
Pauli matrices. After taking the trace and the sum in Matsubara frequencies in @ it is obtained in the zero- 
temperature limit 

where A is an appropriate momentum cutoff to regularize the momentum integral in the ultraviolet, and the quasi- 
particle energy spectrum, e|, is given by 



e e k = ^(e k -efi) 2 + A 2 , e k = V^ + m 2 , e = ±. 



(5) 
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FIG. 1: Energy gap, A, vs the coupling constant g = gA 2 for a free-diquark gas. 
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FIG. 2: Chemical potential, /j,, and mass m vs g = gA 2 for a free-diquark gas. 



The spectra corresponding to different e = ± values denote the particle (e = +) and antiparticle (e = — ) contributions. 

A stable phase must minimize the free energy with respect to the variation of the gap parameter dflo/dA = 0. 
Then, from Q we obtain the gap equation 
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FIG. 3: e^" vs k plotted for different g — gA 2 values. 
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As usual in the study of the BCS-BEC crossover we will consider a canonical ensemble where the particle number 
density, rip = —dClo/dfi, is fixed through the Fermi momentum Pp as np — P f /3tt 2 . Then, from Q we get 



p3 
F 


f d 3 k 




3tt 2 


J A (2tt) 3 


.4 e k. 



with 

Cfe = e fc T M (8) 

Now, we solve numerically the system of Eqs. ((6]) and (O to find the gap A and chemical potential /i, which 
correspond to different values of the coupling constant g. We scale the theory parameters so to guarantee a relativistic 
regime that simulates a quark gas at moderate densities. That is, Pf/A = 0.3, raj A = 0.2. The results for A and 
/i as functions of g, in the interval 0.06 > g > 2, with g = gA 2 , and for A = 602.3 MeV, are shown in Figs. 1 and 2 
respectively. From Fig.l we see that increasing the coupling strength, the energy gap A becomes larger. This suggests 
that the binding energy of the diquark condensate approaches that of a Bose-Einstein condensate (with a smaller 
coherence length £ ~ 1/A) at stronger coupling. To corroborate that this is the case, we should observe the behavior 
of the chemical potential with increasing the coupling strength in Fig. 2. 

As known, the condition /i < m is characteristic of a relativistic Bose gas [19J. From Fig. 2 we see that for this 
simple model there exists a critical value for the coupling constant g cr ~ 1.1 beyond which the condition \i < m is 
satisfied. Hence, the quasiparticle spectrum corresponding to coupling constants smaller and larger than g cr should 
correspond to fcrmion-like and boson-like behaviors, respectively. In Fig. 3, we have plotted the quasiparticle spectra, 
et , corresponding to different values of the coupling constant. The gap, A, and chemical potential, /i, entering in 
the quasiparticle spectrum (JSJ are obtained as solutions of Eqs. © and (|7J| for each value of g. From their graphical 
representations in Fig. 3, we can see that for the spectra corresponding to g — 0.06 and 0.9 the minimum of their 
dispersion relations occurs at k = \J [i 2 — m 2 , with excitation energy given by the gap A, a behavior characteristic 
of quasiparticles in the BCS regime. On the other hand, for g = 2, the minimum of the corresponding spectrum 
occurs at k = 0, with excitation energy — m) 2 + A 2 , which is typical of Bosonic-like quasiparticle. Therefore, it 
is corroborated that g cr is the threshold value for the BCS-BEC crossover in this model. In other words, for g < g cr , 
we have fi > m and the quasiparticles exhibit fermionic-like modes, while for g > g cr , we have fi < m and the 
quasiparticles are characterized by bosonic-like modes. 

We should underline that the obtained value of the critical coupling for the BCS-BEC crossover (g cr ~ 1.1) is in 
the allowed range of values for the strong coupling regime of QCD, where g — gA 2 is expected to be of the order 
of the quark-antiquark coupling G5A 2 [8j. As found in Ref. [18J, once fixed the up and down quark masses with 
equal values m Ut d — 5.5MeV, the four observables of vacuum QCD with values, m T = 135.0 MeV, wik — 497.7 MeV, 
m v > = 957.8 MeV and f„ = 92.4 MeV, are obtained for A = 602.3 MeV (which is the value used in our calculations) 
and G S A 2 = 1.835. 



III. UNSTABLE BEC FREE-DIQUARK REGION 

To find the system EoS it is needed to find the system energy density and pressure. In the case we are investigating 
those magnitudes will depend on the coupling-constant strength. Therefore, varying the values of g from g < g cr to 
g > g cr we will be able to describe the EoS corresponding to the BCS and BEC regimes respectively. 

The energy density and pressure are obtained respectively from the (Too) and (Tu) components of the quantum- 
statistical average of the energy momentum tensor. For an isotropic system, as the one we are considering, the 
covariant structure of the (T^) tensor is given as 20] 

T 

y {T^u) = (^0 + B)g tiv + (fmp + TSju^u,, (9) 

where V is the system volume, T the absolute temperature, S the entropy, and u u the medium 4-velocity with value 
— ^ f 

Up = (1, ) in the rest frame. In ([9]) we introduced the bag constant B to account for the energy difference between 

the perturbative vacuum and the true one. In that way, we are modeling what occurs in the case of quark matter, 

where the asymptotically-free phase of quarks forms a perturbative regime (inside a bag) which is immersed in the 

nonperturbative vacuum. This scenario is what is called the MIT bag model [2l| . The creation of the bag costs free 

energy. Then, in the energy density, the energy difference between the perturbative vaccum and the true one should 

be added. Essentially, that is the bag constant B characterizing a constant energy per unit volume associated to the 

region where the quarks live. From the point of view of the pressure, B can be interpreted as an inward pressure 
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FIG. 4: Energy density, e, and pressure, p, vs the coupling strength g — gA 2 for a free-diquark gas. 



needed to confine the quarks into the bag. In the numerical calculations we will take B 1 / 4 = 145 MeV, which is a 
value compatible with that found in the MIT model. 

Then, the energy density and pressure of the system in the zero-temperature limit are respectively calculated from 

E = VL a + im F + B, p=-n -B (10) 

The results for e and p are plotted in Fig. 4 versus the coupling-constant strength g. There, we can see that the 
system energy density is increasing with the coupling strength, while the pressure is decreasing up to get negative 
values at coupling constants corresponding to the BEC regime. The appearance of a negative pressure for the diquark 
free gas in the BEC region indicates that the free-diquark system is unstable. 



IV. EQUATION OF STATE OF SELF-INTERACTING DIQUARKS 



The system pressure decay in the BEC region obtained in Fig. 4 is an expected result since the absence of repulsion 
between the diquarks makes their Bose-Einstein condensation inevitable with the corresponding decrease of the matter 
pressure. Nevertheless, as we will show in this section, the contribution of the diquark-diquark repulsion in the EoS 
of the strongly interacting system compensates the decreasing tendency due to the Bosc-Einstcin condensation and 
consequently rendering a constant pressure throughout the strongly interacting region. 

The modeling of self-interacting diquarks in the context of a </> 4 boson theory was initially developed in (22j and 
then applied to different situations in (23[. In our case, it can be achieved by introducing a A A 4 term in the free 
energy © 

1 ™ r r1 3 k A 2 
^ = -«E WY Trln [/3G- 1 (^n,k)] + _ + AA 4 , (11) 

Hence, the system energy density and pressure given in (|10j) become 

e = n + AA 4 + /.m F + B, p = -£! - AA 4 - B (12) 

A possible value for the coupling constant A was estimated as A = 27.8 in [22]. It was found taking into account 
the quark interactions in the context of a modified P-matrix formalism of Jaffe and Low [24| . 

The values for A and fi obtained for A = 27.8 from the modified gap equation after including the diquark-diquark 
repulsive interaction term 
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and ||7J), are given in Fig. 5. 

The repulsive interaction between diquarks makes a significant contribution to the energy density and pressure (|12[) 
as can be seen comparing Figs. 4 and 6. From Fig. 6, it is apparent that the instability produced by a negative 
pressure in the BEC region disappears. The matter pressure now remains almost the same in the whole strongly 
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FIG. 5: Energy gap, A, chemical potential n and mass m vs the coupling constant g = gA 2 for a self-interacting diquark gas 
with A = 27.8. 
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FIG. 6: Energy density, e, and pressure, p, vs the coupling strength g — gA 2 for a self-interacting diquark gas with A = 27.8. 



interacting region. In this scenario, the repulsion between diquarks produce enough outward pressure to elude the 
star collapse. On the other hand, this same effect prevents the gas condensation into a zero momentum ground state 
at zero temperature. The fact that the diquark repulsion prevails against the Bose-Einstein condensation is reflected 
in the behavior of the system chemical potential that now can never cross the constant m line in Fig. 5. As we 
checked numerically, once the diquark repulsion is taken into account, the condition fi < m is never reached. Hence, 
the BCS phase is maintained for the whole interaction range. 

If one considers arbitrary values of A, one can show that if the repulsion is weak enough (A < A c = 7.65), the BEC 
dynamics can be reached by increasing g, and consequently, the system pressure becomes negative. On the contrary, 
if A > A c , the pressure is positive for the whole range of considered g-values. But in this last case one has \x > m for 
all those g-values, implying the absence of a BEC region. Notice that A c is smaller than the estimated value A = 27.8 
[2~3 |. In conclusion, we find that there is no way to put together a BEC dynamics with a positive pressure; meaning 
that a gravitational-bound compact star cannot be formed by BEC quark molecules. 

We should highlight that our approach is different from that developed in Ref. 22]. In our case, the diquark 
repulsion effect is treated in the EoS as a dynamical variable, which is determined through the gap equation (|13[) 
and the number density constraint ([7)): while in [22| the contribution of the diquark repulsive potential to the energy 
density and pressure was modeled by assuming an ad hoc Gaussian diquark distribution for the occupation of the 
states with momenta k > 0. In this way, it was prevented that at zero temperature the ground state of the diquark 
gas condensed in the k = state characteristic of the Bose-Einstein condensation phenomenon of a free Boson gas. 



V. SUMMARY AND FINAL REMARKS 



The goal of this paper is to illustrate the behavior of a diquark gas in the strong coupling regime. We started by 
considering a free diquark system. In this system we found that as the strength of the attractive coupling between 
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FIG. 7: Equation of state of the strongly interacting system with diquark-diquark repulsive interaction (A = 27.8) for different 
values of the coupling strength g around the critical value g cr . 



quarks increases, the chemical potential transits from being larger than the quark mass to being smaller, an indication 
of a crossover from the BCS region to the BEC one. A consequence of the crossover to BEC is that the matter pressure 
decays to zero, and even reaches negative values, a sign that the BEC regime cannot be realized in the interior of a 
neutron star. 

We then considered whether this instability could be removed by the introduction of a repulsive force between 
diquarks. In this case, the pressure collapse can be prevented, since a sufficiently strong (A > A c ) diquark-diquark 
repulsion will hinder the overlapping of the diquarks in the ground state. But the implications of reaching a stable 
state in the strong diquark coupling regime is that the system maintains its BCS nature for the whole g- value range. As 
shown then by numerical calculations for A = 27.8 > A c in particular, the contribution of the diquark self-interaction 
is sufficient for stabilizing the system that then acquires an EoS stiff enough (see Figs. 6 and 7) to prevent the collapse 
that would be caused by the formation of a pressureless gas of diquark molecules in the BEC region. The fact that 
the diquark negative pressure prevents the formation of a Bose-Einstein diquark condensate at zero temperature in 
the strongly interacting region is substantiated by the condition jj, > m for the whole interacting region. Thus, a 
self-interacting diquark system will not form a Bose gas. 

The inclusion of the diquark-diquark negative presure in the strongly interacting system makes the EoS stiffcr, 
as can be corroborated from our results. This effect can be important to accommodate quark matter into the EoS 
of neutron stars of high masses, as for instance PSR J1614-2230 [25[ with an inferred value of 1.97 ± 0.04M Q [26| . 
As known, when quark, or other degrees of freedom like hyperons or bosons are considered, the corresponding EoS 
softens and it cannot support highly massive compact stars (27| . unless in the case of quarks, if there exists a color 
superconducting phase with strong interactions between the quarks [281 ] . Nevertheless, the diquark-diquark repulsion 
has not been considered in previous approaches, and its effects are worth to be investigated. 

Nevertheless, we should mention that stars formed by bosons (the so-called boson stars) have been theoretically 
considered since long ago starting with Wheeler's notion of geons [2!| (see Ref. [3(| for recent reviews). A peculiarity 
of those stars is that the internal degrees of freedom are bosonic, so the mechanism to stabilize the star against its 
self-gravity cannot be through the fermionic degeneracy pressure, but by the limitations imposed by the Heisenberg 
uncertainty principle for the mass and radius of the star. 

The result we are reporting in the frame of the simple model of Eq. ([1} should be investigated in more realistic 
models as those of Ref. [1, Q for strong-coupling regimes. Nevertheless, we expect that the main outcome of this 
paper will remain valid. That is, the system pressure will be stabilized by the diquark repulsion. Our expectation 
is based on the fact that by increasing the coupling strength the gap will increase so as to make significant the 
contribution of the outward pressure associated with the diquark repulsive force. A nontrivial problem that remains 
unsolved is to develop the diquark-diquark interaction from first principles. We envision that it will require to start 
from an extended effective theory with a higher number of fermion interactions (as for example, the eight-fermion 
interaction model introduced in Ref. (3l| ) that can give rise in principle to a self- interaction term between the diquark 
condensates. Such a study, however, is out of the scope of the present work. 

An interesting question to be studied in this scenario is the possible effect of an applied magnetic field. As already 
estimated in Ref. 20], magnetic fields of order 10 19 — 10 20 G can coexist in the core of neutron stars with quark 
matter. By increasing the magnetic field strength, A increases [lH, and the system is lead to crossover from the 



BEC region to the BCS one [lj]. At very strong magnetic fields, when all the particles are localized in the lowest 



Landau level, only BCS diquarks are allowed [14| . On the other hand, the pure magnetic contribution to the pressure 
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is negative [33[. Thus, the magnetic field will have a double effect in the pressure whose consequences should be 
elucidated in the frame of the strongly interacting system. 

In those regions of relatively low densities there is of course the possibility that the repulsion between the diquarks 
catalyzes a phase transition to other ground state configurations such as a hadronic phase with a well identified fermion 
nature able to produce the degeneracy pressure needed to compensate for the gravitational pull. Other possibilities to 
be inv estig ated are the viability, through their EoS, of some inhomogeneous phases, such as those formed by density 
waves I34I. quarkyonic chiral spirals foq ]. inhomogeneous Fulde-Ferrel state 36] or quark clusters in solid or liquid 
states [37|, that can in principle be realized in quark matter at moderate density. 
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